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Abstract 

We consider the basic physical properties of matter forming a thin accretion disc in the static 
and spherically symmetric space-time metric of the vacuum f(R) modified gravity models. The 
Lagrangian of the generalized gravity theory is also obtained in a parametric form, and the con- 
ditions of the viability of the model are discussed. The exact Schwarzschild type solution of the 
gravitational field equations in the f{R) gravity contains a linearly increasing term, as well as a 
logarithmic correction, as compared to the standard Schwarzschild solution of general relativity, 
and it depends on four arbitrary integration constants. The energy flux and the emission spectrum 
from the accretion disk around the f{R) gravity black holes are obtained, and they are compared 
to the general relativistic case. Particular signatures can appear in the electromagnetic spectrum, 
thus leading to the possibility of directly testing modified gravity models by using astrophysical 
observations of the emission spectra from accretion disks. 
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I. INTRODUCTION 



Several recent astrophysical observations [If] have provided the astonishing result that 
around 95-96% of the content of the Universe is in the form of dark matter + dark energy, 
with only about 4-5% being represented by baryonic matter. More intriguing, around 
70% of the energy-density may be in the form of what is called "dark energy", and is 
responsible for the acceleration of the distant type la supernovae |2| ■ Hence, today's models 
of astrophysics and cosmology face two severe problems, that can be summarized as the dark 
energy problem, and the dark matter problem, respectively Although in recent years many 
different suggestions have been proposed to overcome these issues, a satisfactory answer has 
yet to be obtained. 

A very promising way to explain these major problems is to assume that at large scales 
the Einstein gravity model of general relativity breaks down, and a more general action 
describes the gravitational field. Theoretical models in which the standard Einstein-Hilbert 
action i8 re p,aced by an arbitrary function of the Ricci scalar R, first proposed in g fiave 
been extensively investigated lately. The only restriction imposed on the function / is to 
be analytic, that is, it must possess a Taylor jseries expansion about any point. Cosmic 
acceleration can be explained 



models have been derived in 



3V f(R) gravity and the conditions of viable cosmo logical 
5J. In the context of the Solar System regime, severe weak 



field constraints seem to rule out most of the models proposed so far 
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although viable 



Uj. The possibility that the galactic dynamic of massive test 



particles can be understood without the need for dar 



framework of f(R) gravity models 
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£ matter was also considered in the 



161 ] . and connections with MOND and 



the Pioneer anomaly further explored by considering an explicit coupling of an arbitrary 

nn 



function of R with the matter 
modified gravity models see 



iagrangian density [171 . Il8j . For a recent review of the f(R) 



The study of the static spherically symmetric vacuum solutions of the gravitational field 
theories is fundamental for the physical understanding and interpretation of the model. In 
particular, the vacuum solutions provide the theoretical basis for the Solar System testing of 
the theories, and for the description of the motion of the test particles around massive bodies. 
It was shown that for a large class of models, including e.g. the f(R) = R — fi A /R model, 
the Schwarzschild-de Sitter metric is an exact solution of the field equations. Solutions in 
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the presence of a perfect fluid were analyzed in [2fJ . Other approaches in searching for exact 
spherically symmetric solutions of f(R) theories of gravity were studied in 21|, respectively. 

Several exact vacuum static and spherically symmetric solutions of the gravitational field 
equations in f(R) gravity were obtained in 23j. The set of the modified Einstein's field 
equations were reduced to a single, third order differential equation, and it was shown how 
one can construct exact solutions in different f(R) models. In particular, a Schwarzschild 
type solution of the field equation was constructed. This solution, containing a term linearly 
increasing with the radial coordinate, as well as a logarithmic term, depends on four new 
arbitrary integration constants, and reduces to the standard general relativistic case by 
imposing the zero value to one integration constant, and by appropriately choosing the 
numerical values of the other constants. However, even that this choice of constants allows 
the model to pass the solar system tests, we expect that at large distances, or in the presence 
of strong gravitational fields, the geometry of the space-time in f(R) modified gravity models 
is different from that of standard general relativity. Therefore, it is important to find a 
method that could allow one to observationally distinguish, and test in an astrophysical 
setting, the possible deviations from Einstein's theory. One such possibility is the study of 
accretion disks around compact objects. 

Most of the astrophysical bodies grow substantially in mass via accretion. Recent ob- 
servations suggest that around almost all of the active galactic nuclei (AGN's), or black 
hole candidates, there exist gas clouds surrounding the central compact object, together 
with an associated accretion disc, on a variety of scales from a tenth of a parsec to a few 
hundred parsecs 25]. These gas clouds are assumed to form a geometrically and optically 
thick torus (or warped disc), which absorbs most of the ultraviolet radiation and the soft 
X-rays. The gas exists in either the molecular or the atomic phase. The most powerful 
evidence for the existence of super massive black holes comes from the VLBI imaging of 
molecular H2O masers in the active galaxy NGC 4258 [26]. This imaging, produced by 
Doppler shift measurements assuming Keplerian motion of the masering source, has allowed 
a quite accurate estimation of the central mass, which has been found to be a 3.6 x 10 7 M Q 
super massive dark object, within 0.13 parsecs. Hence, important astrophysical information 
can be obtained from the observation of the motion of the gas streams in the gravitational 
field of compact objects. 

The mass accretion around rotating black holes was studied in general relativity for the 



3 



first time in [27| . By using an equatorial approximation to the stationary and axisymmetric 
space-time of rotating black holes, steady-state thin disk models were constructed, extending 
the theory of non-relativistic accretion [28|. In these models hydrodynamical equilibrium is 
maintained by efficient cooling mechanisms via radiation transport, and the accreting matter 
has a Keplerian rotation. The radiation emitted by the disk surface was also studied under 
the assumption that black body radiation would emerge from the disk in thermodynamical 
equilibrium. The radiation properties of the thin accretion disks were further analyzed in 



29j and in [30(| , where the effects of the photon capture by the hole on the spin evolution 
were presented as well. In these works the efficiency with which black holes convert rest 
mass into outgoing radiation in the accretion process was also computed. 

Later on, the emissivity properties of the accretion disks were investigated for exotic 
central objects, such as non-rotating or rotating quark, boson or fermion stars 3l|, |32|, 133 ] . 
The radiation power per unit area, the temperature of the disk and the spectrum of the 
emitted radiation were given, and compared with the case of a Schwarzschild black hole of 
an equal mass. 

It is the purpose of the present paper to study the thin accretion disk models applied for 
black holes in f(R) modified gravity models, and carry out an analysis of the properties of 
the radiation emerging from the surface of the disk. 

The present paper is organized as follows. The f(R) gravity generalization of the 
Schwarzschild type solution of general relativity is obtained in Section II. In Section III 
we review the formalism of the thin disk accretion onto compact objects. The basic prop- 
erties of matter forming a thin accretion disc in the space-time metric of the f(R) modified 
gravity models are considered in Section IV. We discuss and conclude our results in Section 
V. In the present paper we use a system of units so that c=G = h = kB = l, where ks is 
Boltzmann's constant. 

II. VACUUM FIELD EQUATIONS IN f(R) GRAVITY 
In the f(R) gravity models the gravitational action is given by 

s = J f(R)V=gd\ (l) 
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where f(R) is an arbitrary functions of the Ricci scalar R. Since we are only interested in 
the vacuum case we do not add a matter Lagrangian to the action. 

Varying the action with respect to the metric g^ u yields the field equations, given by 



F(R)R I1U - ^f(R)g^ - V At V„F(i2) + g^UF(R) = 0, 



(2) 



where we have denoted F(R) = df(R)/dR. Note that the covariant derivative of these field 



22|. 



equations vanishes for all f(R) by means of the generalized Bianchi identities 17, 

In the following we will restrict our study of f(R) gravity models to the static and 
spherically symmetric metric given by 



ds 2 = -e^dt 2 + e x ^dr 2 + r 2 (d9 2 + sin 2 Odcf) 



(3) 



For a metric of the form given by Eq. ([3]), the vacuum field equations of f(R) gravity can 



be expressed as 
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where the prime denotes differentiation with respect to r. Introducing a new variable £ by 
means of the transformation £ = lnr, the field equations Eqs. (j4j) and (jSJ) become 



(i/, € + A >c )F = 0. 



(8) 



1 



+ + \ c)(^« + 2 ) + 2(l-e A ) = - _ A , ., A , 



+ (A e + 4)^, (9) 



where the comma denotes differentiation with respect to the variable £. As a result of 
introducing the new independent variable the basic field equations Eqs. (J4j) and §E§ are 
independent of the coordinate £. It is useful to introduce a formal representation of the 
function F as F t %/F = u, F ^/F = + u 2 , with w a new function of £. Then Eq. (j3J) can 
be written as 



M t + u 



u 



(10) 
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This equation is a Riccati type first order differential equation. By using the function u, Eq. 
(JSJ) becomes 

u K - v c + ^ - -(1/ £ + A, c ) (1/ £ + 2) + 2(1 - e A ) = -2u c - 2u 2 + (A, £ + 4)w. (11) 

Substituting the term u g + w 2 in Eq. ffTTl) with the use of Eq. (jTUJ), we obtain 

2(l-e A ) + z/ iK -z/ iC + ^ + (z/ ie + A if )(l-^ i5 )-(2-z/, c )n = 0. (12) 

The general solution of Eqs. ffTUl) and f[T2"j) gives the general solution of the field equations 
for the static vacuum case of the f(R) gravity models. Once z/(£) and A(£) are specified, 
one can immediately obtain u and then (by integration) F, as well as all the other relevant 
physical quantities. If the function F and the metric tensor coefficients are known, / can 
be obtained as a function of R from Eqs. fl§D and ([7]) in a parametric form, as / = / (£), 
R = R&). 

We now consider general solutions with v £ + A, £ = 0, which we denote as Schwarzschild- 
type solutions. The constant of integration may be set to zero by re-scaling the time coordi- 
nate, so that without a significant loss of generality one may consider the solution u+ A = 0. 
Thus, Eq. (fTUI) reduces to w g + u 2 — u = 0, which provides the solution 

u = ^Tc = 7Tc' (13) 

where C is a constant of integration. Next, from the definition of u, and by reverting back to 
the r coordinate, we obtain F(r) = Ar + B, with A and i? constants of integration. Solving 
Eq. (1121) . one finds the Schwarzschild type general solution in f(R) gravity, given by 
v —A 1 AC 2 C 2 A(B 2 - AC 2 ) 1 f 2 C x \ 2 



where C\ and C 2 are arbitrary constants of integration. An interesting difference to the 
vacuum solutions in general relativity is the presence of the linear term with respect to r, and 
of the term with the logarithmic dependence of r. In the general case the Schwarzschild type 
exact vacuum solution in f(R) gravity depends on four arbitrary integration constants A, B, 
C\ and C 2 , respectively. In order to obtain the Schwarzschild-de Sitter solution of standard 
general relativity, one sets the following values for the constants: A = 0, C 2 = —6BM and 
C\ = 2£> 4 A, respectively, where M is the mass of the central object [23] . 



In the following we assume that G% = —QBM, a condition which is necessary to recover 
the standard Schwarzschild solution. Moreover, we neglect the possible effect of a cosmolog- 
ical constant by taking C\ = 0. By denoting the ratio A/B = 1/r], we represent the vacuum 
metric in f(R) gravity in the form 



, 3M 2M / 6M\ fr\ fry \ / 61, 

e- = e - A = l + 1 + (_)+(_ \ 1 + 1 + ) ln 

r] r V »7 / \Vj \Vj I V V 



V r 



(15) 

By using the Schwarzschild-type solution of the field equations in the modified f(R) 
gravity model given by Eqs. (fT5l) we obtain for the functions f(r), R(r) and dF/dR = 
d 2 f/dR 2 the following expressions 

fir) = —. 1 . - \2Ar] \3 (6M - r) r 2 + rjr (9M + r)) - 3r] 2 (M - r)l 

-6Ar 2 (6M + 77) (r + r?) In [A (l + r/r" 1 )] } , (16) 

i?(r) = — 2 — {77 [6M (2r + 77) (6r 2 + 6r/r - r/ 2 ) + r (-12r 3 - 12r/r 2 + 7r] 2 r + 677 s )] 

rf (r + 77) r 2 

-12r 2 (6M + 77) (r + 77) In [A (l + 777^)] } (17) 

and 

dF_d 2 f_ Ar 3 (r + r,) 3 

dR dR 2 2r 3 + 6r)r 2 + 6t] 2 (r - 2M) 1 ' 

respectively. Generally, the relation between / and R cannot be obtained in an explicit 
analytical form. However, some approximate representations are possible in the limiting 
cases of small and large r, respectively. 

In the limit of small r, so that r is approximately equal to a few M, the logarithmic cor- 
rection In (A + B/r) becomes the dominant term, and we obtain f{r) « — X/ 1 hi-B/r, where 
Xf = r] 2 /6A (6M + 77), and R(r) m ~Xr {t + v) 1 hxB/r, where xr — 12 (6M + 77) /77 s . By 
eliminating r using the relation r = Bexp (—%//) we obtain 

Bexp(-Xff) +V 

If / is small we reobtain the Lagrangian of the standard general relativity, / ~ R. By 
performing a series expansion of the exponential factor we obtain 

Xf BR + xr XfXR V Xr J 
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In the limit of large r, so that r >> M and r >> 77, by keeping only the terms 
in 1/r, we obtain /(r) ~ x/oo + and i?(r) « XRoo/ r ~ 12/r/ 2 , where x/oo = 

-6 [677 + (6M + 77) In A] A and XRoc = 12 [6M77 - r/ 2 - (6M + 77) In A] /r/ 3 . Eliminating 1 /r 
we obtain 

f(R) = — R + A, (21) 

where A = Vl/rfxRoo + x/oo- Therefore at large distances the modified f(R) gravity model 
generates a constant term, the cosmological constant, that is responsible for the accelerated 
expansion of the Universe. 

Eqs. (TT6]) - (TT8]) can be used to discuss the conditions under which the generalized 
Schwarzschild type solution, as well as the corresponding f(R) generalized gravity the- 
ory, represents a viable model. The conditions under which f(R) theories can represent 
viable models of cosmic acceleration have been summarized in [2J]. The requirement of 
the existence of a stable, high curvature regime imposes the constraint d 2 f/dR 2 > 0, for 
R » d 2 f/dR 2 . From Eq. (fl8l) it follows that this condition can be satisfied for all r if 
A < and 77 > 0. Moreover, the condition (r/77) 3 + 3 (r/77) 2 + 3 (r/rj) > 6M/7] must also 
hold for all r, r\ and M. From the tight constraints of the Big Bang Nucleosynthesis and of 
the Cosmic Microwave Background we require that F = df jdR must be a negative, mono- 
tonically increasing function of R that asymptotically approaches to zero from below. The 
condition F = — Ar + B < can be used to constrain the range of the radial variable as 
r > i], where we have assumed that B > 0. 

From the requirement that the effective Newton constant G e // = Gj (1 + F) is not al- 
lowed to change sign we obtain the condition 1 + F > 0, at all finite R. This condition 
gives r < rj + 1/A. Therefore in the present model the range of the radial variable r is 
restricted to 77 < r < rj + 1/A. In order for the model to be viable on very large dis- 
tances, corresponding to cosmic scales, A must have a very small numerical value. Finally, 
F must be a small quantity, a condition which is required for the model to pass the solar 
and galactic scale constraints. In terms of the scalar curvature R, the conditions of viability 
of the f(R) gravity model with Schwarzschild like vacuum solution can be formulated, in 
the small r limit, as —2B (B + rj) /x/Xr < (condition which follows from d 2 f/dR 2 < 0), 
XB./2B <R< xr/ZB+x/Xr/ZB (B + 77), which follows from df/dR < and 1 + df/dR > 0, 
respectively. 

Since A must be a very small (dimensionless) quantity, A << 1, the condition of the 
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smallness of F is automatically satisfied by also assuming a small numerical value for B. By 
assuming that the value of \F\ should not exceed today the numerical value of 10~ 6 at any 
point in the space-time 24|, the requirement that the modified gravity model is viable on 
a scale of r = 1000M, for example, where M is the mass of the central black hole, and by 
representing B as B = j3M, where f3 is a constant, gives the condition (3— 1000 A = 10 _6 /M. 
For a black hole with a mass of around three solar masses we obtain f3— 1000A « 6.74 x 10~ 12 . 
For f3 — 1 it follows that A is very close to A « 10~ 3 . 



III. THIN ACCRETION DISKS ONTO BLACK HOLES 

For the thin accretion disk it is assumed that its horizontal size is negligible as compared 
to its vertical extension, i.e, the disk height H, defined by the maximum half thickness of the 
disk, is always much smaller than the characteristic radius r of the disk, H « r. The thin 
disk is in hydrodynamical equilibrium, where there is only negligible pressure gradient and 
a vertical entropy gradient in the accreting matter. The efficient cooling via the radiation 
over the disk surface prevents the disk from cumulating the heat generated by stresses and 
dynamical friction. In turn this equilibrium causes the disk to stabilize its thin vertical size. 
The thin disk has an inner edge at the marginally stable orbit of the black hole potential, 
and the accreting plasma has a Keplerian motion in higher orbits. 

In steady state accretion disk models, the mass accretion rate Mo is assumed to be a 
constant that does not change with time, and the physical quantities describing the orbiting 
plasma are averaged over a characteristic time scale, e.g. At, over the azimuthal angle 
A0 = 2tt for a total period of the orbits, and over the height H 
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28 



29|. 



The particles moving in Keplerian orbits around the black hole with a rotational velocity 
Q = dip/ dt have a specific energy E and a specific angular momentum L, which, in the steady 
state thin disk model, depend only on the radii of the orbits. These particles, orbiting with 
the four-velocity u M , form a disk of an averaged surface density S, the vertically integrated 
average of the rest mass density po of the plasma. The accreting matter in the disk is 
modelled by an anisotropic fluid source, where the density po > the energy flow vector 
and the stress tensor are measured in the averaged rest-frame (the specific heat was 
neglected). Then the disc structure can be characterized by the surface density of the disk 
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23,1221, 



27|, 



29J. From the 



E(r) = I (po)dz, (22) 

J-H 

with averaged rest mass density (p ) over At and 2it and the torque 

W/= [ (ti)dz, (23) 

J-H 

with the averaged component (t 1 !) over At and 2tt. The time and orbital average of the 
energy flux vector gives the radiation flux T{r) over the disk surface as J-(r) = (g 2 ). 
The stress-energy tensor is decomposed according to 

= p u»u u + 2u^q v) + iT, (24) 

where = 0, u^t^ = 0. The four-vectors of the energy and angular momentum flux are 
defined by —E^ = Tj}(d/dt) v and J M = T^(d/d(p) u , respectively The structure equations 
of the thin disk can be derived by integrating the conservation laws of the rest mass, of 
the energy, and of the angular momentum of the plasma, respectively 
equation of the rest mass conservation, V M (pow M ) = 0, it follows that the time averaged rate 
of the accretion of the rest mass is independent of the disk radius, 

Mo = — 27rrS-u r = constant. (25) 

The conservation law V M -E M = of the energy has the integral form 

[MqE - 2irr{lW/}, r = AnrFE , (26) 

which states that the energy transported by the rest mass flow, MqE, and the energy trans- 
ported by the dynamical stresses in the disk, 2irrQW ( p r , is in balance with the energy radiated 
away from the surface of the disk, AirrFE. The law of the angular momentum conservation, 
V M J M = 0, also states the balance of these three forms of the angular momentum transport, 

[M L - 2mrW/] >r = AirrFL . (27) 

By eliminating from Eqs. (1261) and (12 7p . and applying the universal energy-angular 
momentum relation dE = QdJ for circular geodesic orbits in the form E )T = f2L r , the flux T 
of the radiant energy over the disk can be expressed in terms of the specific energy, angular 



momentum and of the angular velocity of the black hole 27|, |29( , 



T {r) = -Ml „ ^ / ( E _ QL)L rdr . (28) 



A-nr (E - VlL) 2 
10 



In the derivation of the above formula the "no torque" inner boundary condition were 
prescribed, where the torque vanishes at the inner edge of the disk. Thus, we assume that 
the accreting matter at the marginally stable orbit r ms falls freely into the black hole, and 
cannot exert considerable torque on the disk. The latter assumption is valid only if strong 
magnetic fields do not exist in the plunging region, where matter falls into the black hole. 

Another important characteristics of the mass accretion process is the efficiency with 
which the central object converts rest mass into outgoing radiation. This quantity is defined 
as the ratio of the rate of the radiation of energy of photons escaping from the disk surface to 



infinity, and t 



at infinity 



re rate at which mass-energy is transported to the black hole, both measured 
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291 ] . If all the emitted photons can escape to infinity, the efficiency is given 



in terms of the specific energy measured at the marginally stable orbit r ms , 

e = 1 - E ms . (29) 

For Schwarzschild black holes the efficiency e is about 6%, whether the photon capture 
by the black hole is considered, or not. Ignoring the capture of radiation by the hole, e is 
found to be 42% for rapidly rotating black holes, whereas the efficiency is 40% with photon 
capture in the Kerr potential 3p| . 

In order to compute the flux integral given by Eq. (I28p . we determine the radial depen- 
dence of the angular velocity Q, of the specific energy E and of the specific angular momen- 
tum L of particles moving in circular orbits around black holes in the static and spherically 
symmetric geometry given by Eq. ([3]). The geodesic equations for particles orbiting in the 
equatorial plane of the black hole can be written as 

e 2 H 2 = E\ e^ +A V 2 + V eff (r) = E\ r 4 2 = l\ (30) 

where the dot denotes the differentiation with respect to the affine parameter, and the 
effective potential is given by 

Veff(r) =e"(l + ^\ . (31) 
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From the conditions V e ff(r) = and V e ff, r {r) = we obtain 



v_ r e v 



n =\I^Ti ( 32 ) 



E = - , (33) 

~ r 2 Q 

L = - . (34) 

y/e v - r 2 tt 2 

The condition V e ff jrr (r) = gives the marginally stable orbit r ms (or the innermost stable 
circular orbit), which can be determined for any explicit expression of the function u(r). 

It is possible to define a temperature T(r) of the disk, by using the definition of the flux, 
as F{r) = aT 4 (r), where a is the Stefan-Boltzmann constant. Considering that the disk 
emits as a black body, one can use the dependence of T on T to calculate the luminosity 
L (u) of the disk through the expression for the black body spectral distribution 321 ]. 

4 r s rdr 

L (u) = And 2 1 (a;) = - cos^ 3 / — — -, (35) 

7T J r . exp (lo/T) - 1 

where d is the distance to the source, I(u>) is the Planck distribution function, i is the disk 
inclination, and r« and rj indicate the position of the inner and outer edge of the disk, 
respectively. 



IV. THIN DISK ACCRETION PROPERTIES FOR f(R) MODIFIED GRAVITY 
BLACK HOLES 

The effective potential which determines the geodesic motion of the test particles in the 
equatorial plane of the metric given by Eq. (|T5|) can be written as 

O'-O'C-t)-^-?)]}- 

We plot V e ff for different parameters A (or if) in Fig. [TJ Here the parameter A runs over 
the range of —3 x 10 -3 and 3 x 10~ 3 in each plot whereas B is set to M, 2M, 3M and 
AM, so that r) = M/A, 2M/A, 3M/A and AM/ A, respectively. The effective potential of the 
Schwarzschild solution is also plotted for comparison. 
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FIG. 1: The effective potential V e jf(r) for a black hole of total mass M = 2.5 x 10 6 M o and for the 
specific angular momentum L = AM. The metric parameter A is varied in the range of —3 x 10~ 3 
and 3 x 1CT 3 for fixed B in each plot. The parameter B is set to M (the left hand upper plot), 2M 
(the right hand upper plot), 3M (the left hand lower plot) and 4M (the right hand lower plot), 
respectively. The solid line is the effective potential for a Schwarzschild black hole with the total 
mass M. The various values of A are —3 x 1CT 3 (long dashed line), — 1CT 3 (short dashed line), 
10~ 3 (dotted line) and 3 x 1CT 3 (dot-dashed line). 

The unique singular behavior of V e ff (r) at the spatial infinity, as opposed to the asymp- 
totical fall-off of the Schwarzschild potential, is the most striking feature in the upper pairs 
of the plots, where all the curves representing the f(R) potential for different parameters 
tend to — oo as r — > oo. This is only a pure manifestation of the global properties of the 
spacetime for vacuum f(R) gravity models. It is also true for the cases in the lower two 
panels, where the curves plotted above the Schwarzschild limit go to — oo at the spatial 
infinity, after reaching a local maximum at the marginally bound orbit. This indicates the 
obvious fact that though the metric parameters considerably determine the local properties 
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of the gravitational field, the globally hyperbolic geometry must be insensitive to them. 

This behavior of the potential in the f(R) modified gravity models directly follows from 
the presence of the logarithmic correction term in V e ff, and does not depend on the sign of 
the parameters A and B. In the case of a small negative value of A, for some critical value of r 
the argument of the logarithm tends to zero, and the logarithmic term dominates all the other 
terms in the potential, which behaves like V e ff(r) ~ (r/r/) 2 (1 + 6M/77) In (—A + B/r) (the 
minus sign is due to the smallness of the argument of the logarithmic function —A + B/r). 
Interestingly, the same behavior appears if A is a small positive number. In this case for 
B/r — > 0, the potential behaves like V e ff(r) ~ (r/r/) 2 (1 + 6M/77) In A. Due to the smallness 
of A, the potential tends again, in the large r limit, to —00. 

By fixing the parameter B and increasing A from small negative values to zero, we deepen 
the potential well between the marginally stable orbit (indicated by the peak in the left hand 
side), and the marginally bound orbit r m {, (at the local maximum in the right hand side of 
the potential well). At A = we reach the Schwarzschild geometry and for A > we obtain 
deeper potential well, and reduce the local maximum further till it becomes an inflexion 
point. As a result, for any non zero value of A, we can only decrease the radial distance 
between r ms and r m b, as compared to the case of the Schwarzschild potential. Therefore, 
the domain between the two radii for the Keplerian motion of the matter will be decreased 
as well. 

The physical implications of this effect, specific to modified f(R) gravity models, can be 
studied in detail in Figs. [2} where we have plotted the energy flux emitted by the accretion 
disk for the same values of metric parameters A and B. 

The cut-offs of the radial flux profiles in the right hand side at r TO & always appear, no 
matter if the values of A are positive or negative. However, the sign of the values of A 
is already important for the direction of the change of the intensity, as compared to the 
Schwarzschild geometry. Since the potential well is deeper for A > 0, the specific energies of 
the orbiting particles are lower, which decreases the radiated flux over the disk surface. For 
negative values of A we have the opposite effect: in the higher potential well the particles 
have higher specific energies, and more energy is radiated away than in the previous case. 

These effects in the disk radiation can also be observed in the emission spectrum of the 
accretion disk, plotted in Figs. [3j 

In the lower regime, up to the cut-off frequency of the spectrum at about 10 16 Hz, the 
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FIG. 2: The time averaged flux T(r) radiated by the disk for a black hole of total mass M = 
2.5 x 10 6 Mq. The parameter B is set to M (the upper left hand plot), 2M (the upper right hand 
plot), 3M (the lower left hand plot) and 4M (the lower right hand plot), respectively. The solid 
line is the energy flux for a Schwarzschild black hole with total mass M. The various values of 
A are -3 x 1CT 3 (long dashed line), — 10~ 3 (short dashed line), 1CT 3 (dotted line) and 3 x 1CT 3 
(dot-dashed line). 

disk emission is much lower for f(R) black holes than for the Schwarzschild black hole. As A 
tends to zero, the intensity is of course raising up, and approaches the solid line representing 
the case of the Schwarzschild black hole. For frequencies higher than the cut-off value, the 
differences are not so striking, but still observable. For negative values of A, the cut-off in 
the spectrum shifts to higher frequencies, whereas positive values of A produce lower cut-off 
frequency. The effect of varying the parameter B is also important, as all these differences 
are more striking for lower values of B, as can be seen by comparing the upper left hand 
plot, for B = M, with the lower right hand one, corresponding to B = AM. For all the 
cases we have considered here, the smaller values of the emission at lower frequencies have 
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FIG. 3: Emission spectra coL (u) of disks for a black hole of total mass M = 2.5 x 10 6 M o . The 
accretion disk is aligned into a face-on position, with cosi = 1. The parameter B is set to M (the 
left hand upper plot), 2M (the upper right upper plot), 3M (the upper left hand plot) and AM (the 
lower right hand plot), respectively. The solid line is the emission spectrum for a Schwarzschild 
black hole, with the same total mass M. The various values of A are —3 x 1CP 3 (long dashed line), 
— 10~ 3 (short dashed line), 10~ 3 (dotted line) and 3 x 10~ 3 (dot-dashed line). 



the consequence of hardening the disk emission spectrum. The emission spectrum is slightly 
modified due to the effects of varying A, from small negative values to positive ones, which 
softens the emission as we reach positive values of A. 

We also present the conversion efficiency e of the accreting mass into radiation, measured 
at infinity, which is given by Eq. (1291) . for the case where the photon capture by the hole 
is ignored. The value of e measures the efficiency of energy generating mechanism by mass 
accretion. The amount of energy released by matter leaving the marginally stable orbit, 
and falling down the black hole, is the binding energy E ms of the black hole potential. For 
different metric parameters A and B, the values of E ms are given, together with the radii of 
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A 


B [M] 


r ms [M] 


e 







6.0000 


0.0572 


-3 x 10~ 3 


1 


5.8154 


0.0454 


-i(r 3 


1 


5.9218 


0.0531 


10~ 3 


1 


6.1044 


0.0616 


3 x 10~ 3 


1 


6.4706 


0.0712 


-3 x 1(T 3 


2 


5.8861 


0.0510 


-10" 3 


2 


5.978 


0.0551 


10" 3 


2 


6.0489 


0.0593 


3 x 1(T 3 


2 


6.1700 


0.0638 


-3 x 1(T 3 


3 


5.9218 


0.0530 


-10~ 3 


3 


5.9759 


0.0558 


10~ 3 


3 


6.0398 


0.0586 


3 x 1(T 3 


3 


6.1044 


0.0615 


-3 x 10~ 3 


4 


5.9398 


0.0540 


-10~ 3 


4 


5.9850 


0.0561 


1(T 3 


4 


6.0306 


0.0583 


3 x 10~ 3 


4 


6.0766 


0.0604 



TABLE I: The marginally stable orbit and the efficiency for different f(R) black hole geometries. 
The case A = corresponds to the standard Schwarzschild general relativistic black hole. 

the marginally stable orbits, in Table [J where the figures corresponding to the Schwarzschild 
black hole appear in the first line. For A < both the r ms and e are smaller than those for 
the Schwarzschild potential, whereas we obtain higher values than 6 and 0.0572 for A > 0. 
Nevertheless, the efficiency is decreasing with the increasing values of B. For B = 4M the 
values of e are about 10 percent smaller than those for B = M. 
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V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the basic physical properties of matter form- 
ing a thin accretion disc in the Schwarzschild type vacuum space-time metric of the f(R) 
modified gravity models. The physical parameters of the disc - effective potential, flux and 
emission spectrum profiles - have been explicitly obtained for several values of the parameters 
characterizing the vacuum solution of the generalized field equations. All the astrophysi- 
cal quantities, related to the observable properties of the accretion disc, can be obtained 
from the black hole metric. Due to the differences in the space-time structure, the modified 
f{R) gravity black holes present some very important differences with respect to the disc 
properties as compared to the standard general relativistic Schwarzschild case. 

The determination of the accretion rate for an astrophysical object can give a strong 
evidence for the existence of a surface of the object. A model in which Sgr A*, the 3.7 x 
1O 6 M super massive black hole candidate at the Galactic center, may be a compact object 
with a thermally emitting surface was considered in 3^]. For very compact surfaces within 
the photon orbit, the thermal assumption is likely to be a good approximation because 
of the large number of rays that are strongly gravitationally lensed back onto the surface. 
Given the very low quiescent luminosity of Sgr A* in the near-infrared, the existence of a 
hard surface, even in the limit in which the radius approaches the horizon, places a severe 
constraint on the steady mass accretion rate onto the source, M < 10~ 12 M Q yr _1 . This 
limit is well below the minimum accretion rate needed to power the observed submillimeter 
luminosity of Sgr A*, M > 1(T 1O M yr" 1 . 

Thus, from the determination of the accretion rate it follows that Sgr A* does not have 
a surface, that is, it must have an event horizon. Therefore the study of the accretion 
processes by compact objects is a powerful indicator of their physical nature. Since, as one 
can see from Tabled! the conversion efficiency in the case of the f(R) vacuum solutions is 
different as compared to the general relativistic case, the determination of this parameter 
could discriminate, at least in principle, between the different gravity theories, and constrain 
the parameters of the model. 
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